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We have performed numerical simulations of weakly and strongly magnetized relativistic jets em- 
bedded in a weakly and strongly magnetized stationary or mildly relativistic (0.5c) sheath using 
the RAISHIN code. In the numerical simulations a jet with Lorentz factor y = 2.5 is precessed 
to break the initial equilibrium configuration. Results of the numerical simulations are compared 
to theoretical predictions from a normal mode analysis of the linearized RMHD equations de- 
scribing a uniform axially magnetized cylindrical relativistic jet embedded in a uniform axially 
magnetized moving sheath. The prediction of increased stability of a weakly-magnetized sys- 
tem with mildly relativistic sheath flow to Kelvin-Helmholtz instabilities and the stabilization of 
a strongly-magnetized system with mildly relativistic sheath flow is confirmed by the numerical 
simulations. 
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1. Introduction 

Relativistic jets have been observed in galaxies and quasars (AGNs)[|l|, Q], in black hole binary 
star systems (microquasars)[|3|], and are thought responsible for the gamma-ray bursts (GRBs)[U]. 
Proper motions observed in microquasar and AGN jets imply the jet speeds from ~ 0.9c up to 
~ 0.999c. 

While jets at the larger scales may be kinetically dominated and contain relatively weak mag- 
netic fields, the possibility of much stronger magnetic fields certainly exists closer to the acceler- 
ation and collimation region. Recent GRMHD simulations of jet formation performed by Mizuno 
et al. (2006a) [^] indicate that highly collimated high speed jets driven by the magnetic fields 
threading the ergosphere may themselves reside with a broader wind or sheath outflow driven by 
the magnetic fields anchored in the accretion disk. This configuration might additionally be sur- 
rounded by a less collimated accretion disk wind from the hot corona [ph. 

Recent observations of high speed winds in several QSO's (0.1 - 0.4 c) [||] also indicate that 
the highly relativistic jet could reside in a high speed wind or sheath, at least close to the central 
engine. Other observational evidence such as limb brightening has been interpreted as evidence for 
a slower external flow surrounding a faster jet spine A spine-sheath jet structure has also been 
proposed based on theoretical arguments [|T0|]. 

In this paper we investigate the stability properties of highly relativistic jet flows allowing for 
the effects of strong magnetic fields and relativistic flow in a sheath around the highly relativistic jet 
by 3D RMHD numerical simulations. Such 3D numerical simulations of fluid relativistic jets and 
accompanying theoretical work has allowed relatively unambiguous interpretation of the structures 
observed in the numerical simulations] 11, 12, O]. 



2. Numerical Setup 

In order to study the long-term stability of magnetized sheath-spine relativistic jets, we use 
the 3-dimensional GRMHD code "RAISHIN" with Cartesian coordinates in special relativity. The 
method is based on a 3+1 formalism of the general relativistic conservation laws of particle number 
and energy momentum, Maxwell equations, and Ohm's law with no electrical resistance (ideal 



MHD condition) on a curved spacetime [|14|]. The RAISHIN code can perform special relativistic 
computations in Minkowski spacetime by changing the metric. 

We consider the following initial conditions for the simulations: "preexisting" jet flow is estab- 
lished across the computational domain. This represents the case in which a leading Mach disk and 
bow shock has passed. The jet flow is surrounded by a low-density external wind (medium). For all 
simulations, the ratio of densities is pj/p e = 2.0, where p is the mass density in the proper frame. 
The jet flow has Vj = 0.9165c and y = (1 — v 2 ) -1 / 2 = 2.5. The initial magnetic field is assumed to 
be uniform and parallel to the jet flow. We have performed two sets of simulations. In the weakly 
magnetized (RHD) simulations, the relevant sound speeds are a e = 0.574c and aj = 0.511c, the 
relevant Alfven speeds are VAe = 0.0682c and vaj = 0.064c, and the Alfven speed is much smaller 
than sound speed. In the strongly magnetized (RMHD) simulations. The relevant sound speeds 
are a e = 0.30c and aj = 0.226c, the relevant Alfven speeds are VAe = 0.56c and vaj = 0.45c, and 
the Alfven speed is comparable to the sound speed. In order to investigate the effect of an external 
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Figure 1: Three-dimensional isovolume image (a) and transverse (b) cut (z = 30Rj) of the weakly magne- 
tized case without an external wind. The color scales show the logarithm of density. The white lines indicate 
magnetic field lines. The arrows depict transverse velocities. 

wind, we have performed simulations with no external wind (v e = 0.0c) and a mildly relativistic 
external wind (v e = 0.5c) for both weakly- and strongly-magnetized cases. The computational do- 
main is 3Rj x 3Rj x 60Rj with 60 x 60 x 600 computational zones (10 computational zones span 
Rj). We impose outflow boundary conditions on all surfaces except the inflow plane, z = 0. A 
precessional perturbation of angular frequency coRj/vj = 0.93 is applied at the inflow by imposing 
a transverse component of velocity with v± = O.Olvy. The simulations are halted after ~ 60 light 
crossing times of the jet radius,and before the perturbation has crossed the entire computational 
region. 

3. Results 

Figure 1 shows the structure of a weakly magnetized jet (three-dimensional isovolume image 
and transverse cuts) without an external wind at t = 60. The precession at the jet inflow plane 
excites the Kelvin-Helmholtz (KH) instability. The initial perturbation propagate down the jet and 
grows as a helical structure. Beyond Rj = 40 the KH instability disrupst the jet structure strongly. 
The magnetic field is strongly distorted and becomes complicated. The transverse cut at z = 30Rj 
seen Fig lb, shows strong interaction with the external medium at this distance. Transverse veloc- 
ities show circular motions near the jet surface caused by the helically twisted jet. 

To investigate simulation results quantitatively, we take one-dimensional cuts through the 
computational box parallel to the z-axis at different radial distances along the transverse x-axis 
at x/Rj = 0.2, 0.5, and 0.8, as shown in Figure 2. 

The weakly magnetized simulation results show velocity oscillation from the growing helical 
KH instability. The simulations are fully evolved to a non-linear phase by about 32Rj for no exter- 
nal wind and beyond about 3SRj for an external wind. The external wind has reduced the growth 
of KH instability and has delayed the onset of the non-linear phase. The dominant wavelength 
of oscillation is X /Rj ~ 6 in both weakly magnetized cases. The strongly magnetized simulation 
results without an external wind also show growing velocity oscillation (Figure 2c), albeit more 
slowly. The simulation is still evolving linearly and has not reached the non-linear phase. From 
the comparison with the weakly-magnetized cases, the stronger magnetic field reduces the growth 
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Figure 2: Radial velocity (v A ) along the one-dimensional cuts parallel to the jet axis located at x/Rj = 0.2 
(solid line), 0.5 (dotted line) and 0.8 (dashed line) for the weakly (a, b) and strongly (c, d) magnetized cases 
without an external wind (a, c) and with an external wind (b, d). 



rate of the helical KH instability. The strongly magnetized simulation result with an external wind 
reveals a damped oscillation. Therefore the external wind in the strongly magnetized case leads to 
damping of KH instability and stabilizes the jet. The dominant oscillation wavelength is k/Rj ~ 5 
in both of the strongly magnetized cases. 

The simulations confirm theoretical predictions for the stability of a jet spine-sheath or jet- 
wind configuration Hardee (2006) [15]. Theoretical predictions arise from a normal mode analysis 
of the linearized ideal RMHD equations describing a uniform axially magnetized cylindrical rel- 
ativistic jet embedded in a uniform axially magnetized moving sheath. The most important theo- 
retical prediction obtained from the normal mode analysis is that the presence of an external wind 
and strong magnetic field can reduce the growth rate or lead to damping of the KH instability. This 
prediction is confirmed by the present simulation results. 



4. Conclusion 



We have performed numerical simulations of weakly and strongly magnetized relativistic jets 
embedded in a weakly and strongly magnetized stationary or mildly relativistic (0.5c) sheath using 
the RAISHIN code 014}]. In the numerical simulations a jet with Lorentz factor 7 = 2.5 is precessed 
to break the initial equilibrium configuration. Results of the numerical simulations have been 
compared to theoretical predictions from a normal mode analysis of the linearized RMHD equa- 
tions describing a uniform axially magnetized cylindrical relativistic jet embedded in a uniform 
axially magnetized moving sheath. The prediction of increased stability of weakly-magnetized 
systems with mildly relativistic sheath flow to Kelvin-Helmholtz instability and the stabilization 
of a strongly-magnetized system with mildly relativistic sheath flow is verified by the numerical 
simulation results. 
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